Relativistic acceleration of Landau resonant particles as a consequence of 
Hopf bifurcations. 
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Using bifurcation theory on a dynamical system simulating the interaction of a particle with an obliquely 
propagating wave in relativistic regimes, we demonstrate that uniform acceleration arises as a consequence of 
Hopf bifurcations of Landau resonant particles. The acceleration process arises as a form of surfatron estab- 
lished through the locking in pitch angle, gyrophase, and physical trapping along the wave-vector direction. 
Integrating the dynamical system for large amplitudes {SB/Bq ^ 0.1) obliquely propagating waves, we find 
that electrons with initial energies in the keV range can be accelerated to MeV energies on timescales of 
the order of milliseconds. The Hopf condition of Landau resonant particles could underlie some of the most 
efficient energization of particles in space and astrophysical plasmas. 
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I. INTRODUCTION 

The main purpose of this Letter is to study the 
wave-particle interaction in the general case of oblique 
propagation and in the relativistic limit. The more 
recent inclusion of oblique waves in modeling is not 
fortuitous^ Obliquely propagating waves are not 
only theoretically predictecP but are also observed in 
numerous and diverse space plasma regimes, e.g., in 
fast streams of the solar wincpl^ upstream of the bow 
shock^, magnetosheatlF', cometary environmeniP, the 
Earth's plasma sheet boundarjEl and the Van Allen 
radiation bell!^. The interest in obliquely propagating 
waves mainly resides in the inclusion of an electric field 
along the background magnetic field. Unlike the parallel 
propagation case and due to constraints imposed by 
Maxwell-Lorentz invariant quantities, it is impossible 
to find an inertial frame of reference for which the 
electric field vanishe^^i'. The electric field of the wave 
can therefore provide acceleration mechanisms and/or 
physical trapping along the background magnetic field. 
Similarly, the study of relativistic regimes to modeling 
space and astrophysical plasmas now appears more than 
ever as absolute necessity, and more so for the numerous 
problems (i.e. cosmic rays, radiation belt electrons) 
where the acceleration or/and injection of charged 
particles to relativistic energies remains misunderstood. 
In the following report, we ignore the more commonly 
studied case of cyclotron-resonant particles (cj ^ Q) 
and concentrate instead on the parameter space where 
the Landau resonance (w ^ 51) is accessible. Even 
though cyclotron-resonant interactions are observed and 
understood to play an important role in the kinetic 
description of space and astrophysical plasmas, other 
means of energy exchange between waves and parti- 
cles have been lurked behind a predominant focus on 



cyclotron-resonance. For instance, the perpendicular 
heating in fast streams of the solar wind could as well re- 
sult from trapped particles caught in broadened Landau 
resonance ins tead of the commonly assumed cyclotron 
resonanc d^ ^ ^^ l This report also aims at providing further 
ground for such views, but in the context of relativistic 
and weakly collisional plasmas. 



II. DYNAMICAL SYSTEM 

This problem is addressed by using dynamical systems' 
theory, which although lacking levels of self-consistency 
that simulations provide, can facilitate the understand- 
ing of complex systems such as plasmas and provide for 
an intuitive leap between theoretical models and simula- 
tions. Hence, our study begins by writing the equation of 
motion for a particle in an electromagnetic field as follow 



dp 



E(x,i) 



X B(x,t) 



(1) 



for a particle of momentum p = rwyv, rest mass 
TO, charge e and Lorentz contraction factor 7 — 
\/l + p^/m^c^. The fields topology consist of an elec- 
tromagnetic obliquely propagating wave of amplitude 
((5E, (5B) superposed on a background magnetic field Bq. 
We choose the electromagnetic wave vector k to point 
in the z direction and the background magnetic field to 
lie in the y — z plane. Hence, the propagation angle, 
9, denoting the obliqueness of the wave, is defined as 
k • Bq = kBo cos{9). The magnetic field components of 
the wave are written as 



SB^ = 6B sm{kz — ujt) 
SBy = SBcos{kz — Lot), 



(2) 



''J Electronic mail: 



a.osmane@unb.ca 



with the electric components provided by Faraday's law, 
ck X (5E(k, w) = aj(5B(k, w). We can then write the 
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dynamical system equations for the chosen electromag- 
netic field topology in terms of the following variables : 
v<j, ~ Lo/k, ~ mjv^, Q.I = eSB/mcj, Uq — eBQ/mcj. 
Hence, obtaining the following set of equations : 

' Px — Py^o cos{0) + (p$ — Pz)^i cos(fcz — Lut) + Pz^o sin(6') 
Py = —px^o cos(0) + {pz — sin(fc2; — ujt) 

Pz = —px^Q sin(0) + Px^i cos(fcz — ujt) — Pyf^i sin(fcz — ujt) ^ 

Z = PzV^/p^ 

(3) 

It follows that the dynamical gyrofrequencies (r2o,i^i) 
can be tracked as the fifth variable of the dynamical sys- 
tem : 



50 



■ d feBo 
at \mcj 



-no 



pc 



rn^c^ + p'^c^ 



P- 



(4) 



We now proceed by eliminating the time dependence by 
making the following change of variables : 



Px ^Px,Py ^Py,Pz 



{Pz~P4,), z' ^{z~v^t). (5) 



Hence, we write the the dynamical system in terms of 
the primed variables as follow : 

p'x = ^opy cos{9) - Vlip'^ cos{kz') + f^o(Pz +Pct>) sin(6') 
p'y = — iloP^ cos{9) + flip'^ sin(fcz') 

p', = -ilop'x sm{9) + n,{^){p'^ cos(fcz') - p'y sin(fcz')) 

Z' = p'zV^/p^ 

(6) 

with the refractive index n? — (? Iv\. In addition to 
equation (4), one can compute the particle orbits for a 
class of parameters 9,n,Si = r^o/f^i and 62 = w/fioT- 



III. HOPF BIFURCATIONS OF LANDAU RESONANT 
TRAPPED ORBITS— 

Despite the apparent simplicity of the dynamical sys- 
tem, the inclusion of the relativistic terms results in a 
number of interesting properties that extend beyond the 
scope of this report. We hereafter focus on one of these 
properties arising from the bifurcation in stability of fixed 
(stationary) points. Indeed, it can easily be shown that 
the set of equations (6) possesses a class of fixed (station- 
ary) points that can be represented as follow : 



70 



1 



Pxo ^ PzQ — 0; 
1 

f (l + tan2(0))' 



p'yQ = --P* tan(6') 
Z^kz' = 0,7r. 



These points are informative of the values for which 
a particle is physically trapped by the electromagnetic 
field. Their values in velocity space correspond to the 
Landau resonance condition. This can be more clearly 
seen if one applies the inverse of the translation in equa- 
tion (5) followed by a rotation in a system of coordi- 
nate with the z axis parallel to the background mag- 
netic fielcf^. The next fundamental step in dynamical 
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FIG. 1. Eigenvalues' dependence on the propagation angle d 
for fixed parameters 5i =0.1, &2 = 0.0696, = 2 and the 
fixed point of component Zo = 0. The bifurcation through the 
positive real axis takes place for propagation angle Qc. = 60°. 
The fixed point is stable for 9 < 9^ and unstable for 9 > 9c- 



system theory is to investigate the stability of the fixed 
points. In order to do so, we apply a basic Lyapunov 
linear analysis that can be found in any textbook on dy- 
namical systemJi^l Hence, solving the eigenvalue prob- 
lem (J — AI) = for the Jacobian J and eigenvalue A, we 
find a bi-quadratic polynomial function in A that can be 
written as x(A) = A'^ + r^iA^ + 772 = 0, with the constant 
coefficients 771 and 772 given by the following expressions 



5i n 
(5270 n 
Si 



tan(e') 



cos2(6') 



V2 



61 



Sill 



±2sin(6l) 



■ sin(6') cos(0). 



s\n^{9) sin(6l) 
5i 'n? 



(7) 



with the ± symbol denoting the values for Z = and 
Z = TT. A close look at the coefficients of equation set 
(7) shows that all four eigenvalues will cross the zero real 
axis when the condition 



- 1 = tan^ (61) 



(8) 



is respected. That is, for parameter values corresponding 
to 7q"^ =0 and resulting in A^ ~ 0. This condition can 
be more clearly expressed through Figure 1 where we 
plotted the dependence of the real and imaginary part 
of all four eigenvalues as a function of 9, while keeping 
the remaining parameters ((5i,(52,n^) constant. We ob- 
serve that when the condition in equation(8) is respected, 
the equilibrium evolves from stable to unstable equilib- 
rium since the real part of one of the eigenvalues becomes 
positive. This type of bifurcation, where pairs of com- 
plex conjugate eigenvalues cross through the imaginary 
axis, is the well-known Hopf bifurcation^^. The fixed 
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FIG. 2. Particle orbits for parameters Si — 0.1, S2 = 0.0696, 
— 4:, 6 — 6c — 1° and initial conditions w^q = 0, v'yQ = 
— Ucj tan(S) — 1.6u$, v'^q = —v<f, Zq = 0. 




FIG. 3. Particle orbits for parameters 5i — 0.1, 82 = 0.0696, 
= 4:, 9 — 6c. The orbit is locked in phase-space and trapped 
along Z. 



point for Z — n is linearly unstable for every parameter 
range and values chosen in Figure 1 with the exception 
of the parameters for which equation (8) is respected. 
The eigenvalue profile and the following conclusions do 
not differ significantly for low-frequencies (^2 < 1) and 
large-amplitudes {Si < 1), that is the relevant range of 
parameters in space and astrophysical plasmas. 



IV. NUMERICAL INTEGRATION— 

We now investigate the effects of the bifurcation on the 
particles belonging to the basin of attraction of the fixed 
point with component Z = 0. For the sake of clarity 
and in order to underlie the physical processes at play, 
we choose to represent the particle momentum in terms 




FIG. 4. Particle orbits for parameters 5i =0.1, ^2 ~ 0.0696, 
ri^ = i, 6 = dc- The orbit is uniformly accelerated once locked 
in phase-space. 




FIG. 5. Particle orbits for 6 > dc- A torus emerge from the 
two Hopf bifurcations. 

of spherical coordinates (p', a, $) instead of the carte- 
sian coordinates {p'^,p'y,p'^) in equation set (6). The 
transformation from one expression of the momentum 
through the other can be made using the following def- 
initions for the magnitude p' = y^P^^+KF~l-pl^, the 
pitch angle tan(Q!) — and the dynamical gyrophase 

tan($) = ^ - p;eosW+p',sin(9) ^ ^^i^re the parallel and 
perpendicular symbols are with respect to Bq. In Figure 
2 the orbit of a particle interacting with a large ampli- 
tude and low- frequency wave is shown for a, $ and Z 
for a propagation angle nearly obeying the condition of 
equation (8). For such parameter values, the fixed point 
is linearly stable and a particle belonging to the basin 
of attraction will remain physically trapped along the 
wave-vector. The orbit eventually closes onto itself as 
the particle bounces back and forth in the wave poten- 
tial, alternatively loosing and gaining energy with no net 
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FIG. 6. Lorentz factor 7 as a function of time for Si = 
(0.04,0.05,0.06), S2 = 0.1, n'^ ^ 9, 6 = Oc and initial condi- 
tion — 0, Vy = —11$ tan(Sc) — fs/S, Z = fca' = 0. The 
particle can be accelerated to MeV energies for time scales of 
less than a millisecond. 



gain over one period. If we modify the propagation angle 
such that the condition in equation (8) is respected, we 
can see from Figures 3 and 4 that a particle will asymp- 
totically converge into a point in the (a, $) phase-space 
while it gets trapped along Z and subsequently diverges 
to infinity in momentum. That is, the fixed point be- 
comes an attractor along (a,$), and the particles ini- 
tially belonging to the basin of attraction will be locked 
forever in phase-space and experience uniform accelera- 
tion through the constant electric field observed by the 
particle. If we increase the propagation angle further, 
such that > ^c, we find that two dimensional torus are 
created (see Figure 5) and that particles can be neither 
physically trapped nor uniformly accelerated. However, 
even though the dynamics for 9 > Oc deserves a study 
of its own, we now focus on the case 9 ^ 9c for which 
particles can be energized irreversibly. 



V. DISCUSSION— 

Phase-locking and trapping of a particle while submit- 
ted to a constant electric field has been predicted and 
observed before in previous models using Hamiltonian 
and/or asymptotic approaches. Indeed, similar processes 
for a wide range of electrostatic and/or electromagnetic 
topologies have been referred to as surfatron^*'^^ . The 
novel result reported here is that the mechanism arises 
as a result of Hopf bifurcations at the Landau resonance. 
Qualitatively, one could describe this result by the state- 
ment that there are resonances of Landau resonances re- 
sulting in efficient energization. 

In order for this mechanism to be reproducible in space 
and astrophysical plasmas, the acceleration should take 
place for realistic wave amplitudes, and belong to a suffi- 



ciently wide volume of velocity-space to affect a portion 
of a distribution functions. In other words, the basin of 
attraction should be wide enough, and particles of moder- 
ate energies should be able to be attracted into it. Choos- 
ing a highly oblique wave 9 ~ 71°, with a wave frequency 
given by S2 — O-l and amplitudes of the order of a per- 
cent of the background magnetic field, we integrate the 
dynamical system for a few wave periods. The results are 
shown in Figure 6 for different wave amplitudes. Whereas 
the irreversible acceleration, coinciding with the Hopf bi- 
furcations in phase space, requires a wave amplitude of 
61 > 0.06 capable of physical trapping along Z; particles 
can still be brought to relativistic energies if they are 
caught in sufficiently close to the fixed point. The size of 
the basin of attraction varies significantly as a function of 
the parameter S2 but we nonetheless find that particles 
with moderate energies ranging from few keV to hun- 
dreds of keV, can be accelerated to relativistic levels on 
timescales comparable to the gyroperiod. The basin of 
attraction is therefore wide enough to affect a significant 
portion of a distribution function. In the case shown in 
Figure 6 for ^2 = 0.1, electrons with few hundreds keV 
are accelerated to MeV energies on timescales less than 
a millisecond. Such timescales of energization for large- 
amplitudes Si ^ 0.06 and low frequency 62 ^ 0.1 sug- 
gest that this mechanism could be of interest for studies 
of space (e.g. radiation belts) and cosmic (e.g. galactic 
plasma) plasmas if permeated by large-amplitude oblique 
waves. Further work is now currently underway to apply 
methods presented in this paper to the specific problem 
of electron acceleration in the planetary radiation belts. 
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